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ABSTRACT 

We fit every emission line in the high-resolution Chandra grating spectrum of 
C Pup with an empirical line profile model that accounts for the effects of Doppler 
broadening and attenuation by the bulk wind. For each of sixteen lines or line com- 
plexes that can be reliably measured, we determine a best-fitting fiducial optical depth, 
t* = kM/47tR >i ,v 00 , and place confidence limits on this parameter. These sixteen lines 
include seven that have not previously been reported on in the literature. The ex- 
tended wavelength range of these lines allows us to infer, for the first time, a clear 
increase in r* with line wavelength, as expected from the wavelength increase of bound- 
free absorption opacity. The small overall values of t* , reflected in the rather modest 
asymmetry in the line profiles, can moreover all be fit simultaneously by simply as- 
suming a moderate mass-loss rate of 3.5 ± 0.3 x 10~ 6 Mq yr" 1 , without any need 
to invoke porosity effects in the wind. The quoted uncertainty is statistical, but the 
largest source of uncertainty in the derived mass-loss rate is due to the uncertainty 
in the elemental abundances of £ Pup, which affects the continuum opacity of the 
wind, and which we estimate to be a factor of two. Even so, the mass-loss rate we find 
is significantly below the most recent smooth-wind Ha mass-loss rate determinations 
for £ Pup, but is in line with newer determinations that account for small-scale wind 
clumping. If £ Pup is representative of other massive stars, these results will have 
important implications for stellar and galactic evolution. 

Key words: stars: early-type - stars: mass-loss - stars: winds, outflows - stars: 
individual: £ Pup - X-rays: stars 



1 INTRODUCTION 

Massive stars can lose a significant fraction of their orig- 
inal mass during their short lifetimes due to their strong, 
radiation-driven stellar winds. Accurate determinations of 
these stars' mass-loss rates are therefore important from an 
evolutionary point of view, as well as for understanding the 
radiative driving process itself. Massive star winds are also 
an important source of energy, momentum, and (chemically 
enriched) matter deposition into the interstellar medium, 
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making accurate mass-loss rate determinations important 
from a galactic perspective. 

A consensus appeared to be reached by the late 
1990s that the mass-loss rates of O stars were accurately 
known observationally and th eoretically, using th e modified 
(|Pauldrach et al.lll986t ) CAK (|Castor et al.lll975l 1 theory of 
line-driven stellar winds. This understanding was thought to 
be good enough that UV observations of spectral signatures 
of their winds could be used to determine their luminosi- 
ties with sufficien t accuracy to m ake extragalactic O stars 
standard candles l|Puls et alj|l99fih ■ 

This consensus has unraveled in the last few years, 
mostly from the observational side, where a growing ap- 
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preciation of wind clumping - a n effect whose i mpor- 
tance has long been recognized dEversberg et all 1 19981 ; 
iHillier fe Miilerfil999l : lHamann fe Koesterkdl 19991 ) - has led 
to a re-evaluation of mass-loss rate diagnostics, including Ha 
emiss ion, radio and IR f r ee-free emission, and UV absorp- 
tion l)Bouret et al ] 120051 : IPuls et al ] 120061 : iFullerton et al] 
2006). Accounting for small-scale clumping that affects den- 
sity squared emission diagnostics - and also ionization bal- 
ance and thus ionic column density diagnostics like UV reso- 
nance lines - leads to a downward revision of mass-loss rates 
by a factor of seve ral, with a fair amount of controversy over 
the actual factor llHamann et al.ll200Sl : [Puis et al]l20Qgh . 

X-ray emission line profile analysis provides a good and 
independent way to measure the mass-loss rates of O stars. 
Like the UV absorption line diagnostics, X-ray emission pro- 
file diagnostics are sensitive to the wind column density 
and thus are not directly affected by clumping in the way 
density-squared diagnostics are. Unlike the UV absorption 
line diagnostics, however, X-ray profile analysis is not very 
sensitive to the ionization balance; moreover, as it relies on 
continuum opacity rather than line opacity, it is not sub- 
ject to the uncertainty associated with saturated absorption 
lines that hamper the interpretation of the UV diagnostics. 

In this paper, we apply a quantitative line profile anal- 
ysis to the Chandra grating spectrum of the early O super- 
giant, £ Pup, one of the nearest O stars to the Earth and 
a star that has long been used as a canonical example of 
an early O star with a strong radiation-driven wind. Previ- 
ous analysis of the same Chandra data has established that 
the kinematics of the X-ray emitting plasma, as diagnosed 
by the line widths, are in good agreement with wind-shock 
theory, and that there are modest signatures of attenua- 
tion of the X-rays by the dominant cold wind component in 
which the shock-hea ted X-ray emitting plasma is embedded 
jKramer et alj|2003h . 

The work presented here goes beyond the profile anal- 
ysis reported in that paper in several respects. We analyze 
many lines left out of the original study that are weak, but 
which carry a significant amount of information. We better 
account for line blends and are careful to exclude those lines 
where blending cannot be adequately modelled. We model 
the continuum emission underlying each line separately from 
the line itself. We use a realistic model of the spectrometers' 
responses and the telescope and detector effective area. And 
we include the High Energy Grating (HEG) spectral data, 
where appropriate, to augment the higher signal-to-noise 
Mediu m Energy Grating (MEG) data that iKramer et al] 
(2003) reported on. 

Implementing all of these improvements enables us to 
derive highly reliable values of the fiducial wind optical 
depth parameter, t* = «M/47rR*Voo, for each of sixteen 
emission lines or line complexes in the Chandra grating spec- 
trum of ( Pup. Using a model of the wavelength-dependent 
wind opacity, «, and values for the star's radius, R*, and 
wind terminal velocity, Voo, derived from UV and optical 
observations, we can fit a value of the mass-loss rate, M, to 
the ensemble of t* values, and thereby determine the mass- 
loss rate of f Pup based on the observed X-ray emission line 
profiles. 

In doing this, we also can verify that the wavelength- 
dependence of the optical depth values - derived separately 
for each individual line - is consistent with that of the 



atomic opacity of the bulk wind, rather than the gray opac- 
ity that woul d, for example, be obtained from an extremel y 
porous wind ijOskinova et al.ll2006l : lOwocki fe CoherJl200d ). 
While a moderate porosity might reduce somewhat the 
effective absorption while still retaining some wavelength 
dependence, for simplicity our analysis here assumes a 
purely atomic opacity set by photoelectric absorption, with 
no reduction from porosity. This assumption is justified 
by the large porosity lengths required for any apprecia- 
ble p orosity effect on line profile shapes (|Owocki fc Cohen] 
2006) and the very small-scale clumping in state-of-the- 
art two-dimensional rad iation hydrodynamics simulations 
IIDessart fc Owockj|2003h . Furthermore, preliminary results 
indicate that profile models that explici tly include porosit y 
are not favored over ones that do not (|Cohen et alj fefJOS). 
We will extend this result in a forthcoming paper but do 
not address the effect of porosity on individual line profile 
shapes directly in the current work. 

The paper is organized as follows: We begin by describ- 
ing the Chandra data set and defining a sample of well be- 
haved emission lines for our analysis in §2. We briefly eval- 
uate the stellar and wind properties of £ Pup in §3. In §4 
we describe the empirical profile model for X-ray emission 
lines and report on the fits to the sixteen usable lines and 
line complexes in the spectrum. We discuss the implications 
of the profile model fitting results in §5, and summarize our 
conclusions in §6. 



2 THE Chandra GRATING SPECTRUM 

All the data we use in this paper were taken on 28-29 
March 2000 in a single, 68 ks observation using the Chandra 
High-Energy Transmission Grating Spectrometer (HETGS) 
in conjunction with the Advanced CCD Im aging Spectrom- 
eter ( ACIS) detector in spectroscopy mode (jCanizares et al.l 
2005). This is a photon counting instrument with an ex- 
tremely low background and high spatial resolution (~ 1"). 
The first-order grating spectra we analyzed have a total of 
21,684 counts, the vast majority of which are in emission 
lines, as can be seen in Fig. [T] We modelled every line or 
line complex - 21 in total - as we describe in §4, and indi- 
cate in this figure which of the lines we deemed to be reliable. 
We only include lines in our analysis that are not so weak 
or severely blended that interesting parameters of the line- 
profile model cannot be reliably constrained. (See ^4.2.3l for 
a discussion of the excluded line blends.) 

The HETGS assembly has two grating arrays - the 
Medium Energy Grating (MEG) and the High Energy Grat- 
ing (HEG) - with full-width half maximum (FWHM) spec- 
tral resolutions of 0.0023 A and 0.0012 A, respectively. This 
corresponds to a resolving power of 1Z ~ 1000, or a velocity 
of 300 km s" 1 , at the longer wavelength end of each grating. 
The wind-broadened X-ray lines of £ Pup are observed to 
have Wfwhm ~ 2000 km s _1 , and so are very well resolved 
by Chandra. The wave length calibration of the HETGS is 
accurate to 50 km s" 1 [(Marshall et al.ll2004h . 

The two gratings, detector, and telescope assembly have 
significant response from roughly 2 A to 30 A, with typical 
effective areas of tens of cm 2 , which are a strong function 
of wavelength. In practice, the shortest wavelength line with 
significant flux in the relatively soft X-ray spectra of O stars 
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Figure 1. The entire usable portions of the MEG (top) and HEG (bottom) first order (negative and positive orders coadded) spectra 
of £ Pup. The binning is native (5 mA for the MEG and 2.5 mA for the HEG). Vertical dashed lines in the data panels themselves 
represent the laboratory rest wavelengths of all the 21 lines and line complexes we fit with the profile model. Solid (red) vertical lines 
between the two spectral plots indicate the lines we successfully fit with profile models and lines we attempted to fit but which were too 
blended to extract meaningful model parameters are indicated by dashed (green) lines. For all blended emission lines we show only one 
of these solid or dashed lines between the panels, and align it with the bluest line in the blend. 



like £ Pup is the S xv line complex near 5 A, and the longest 
wavelength feature is the N VII Lye* and N VI He/3 line blend 
at 24.781, 24.890 A. The HEG response is negligible for lines 
with wavelengths longer than about 16 A. 

The X-ray spectrum of £ Pup consists of emission lines 
from H-like and He-like ionization stages of N, O, Ne, Mg, 
Si, and S, and numerous L-shell lines of iron, primarily Fe 
xvil. The Lya lines and often the /? and even 7 lines of 
the Lyman series are seen for the H-like ions. There is a 
weak bremsstrahlung continuum beneath these lines. Over- 
all, the spectrum is consistent with an optically thin, thermal 
plasma in ionization equilibrium with a range of tempera- 
tures from one to several million degrees present. It is pos- 
sible that there are deviations from equilibrium, although 
the spectrum is not of high enough quality to show this. 
There is some evidence from the XMM-Newton RGS spec- 
trum that a few of the emission lines are optically thick 
l|Leutenegger et al] [2007); a possibility we will take into ac- 
count when discussing the results for those lines. 



3 THE STAR 

£ Puppis is a relatively nearby, bright, and well-studied early 
O supergiant (04 If) that shows the enhanced nitrogen and 
deficient carbon and oxygen that is indicative of CNO cycle 
proce ssed material. Helium is also overabundant ()Puls et al.l 
[2006J). The star's rapid rotation may explain the photo- 
spheric abundances, though they may instead have resulted 
from the supernova explosion that is invoke d to explain its 
high space velocity (|Vanbeveren et all 1 1 9981 '). On the other 
hand, no special mechanism may need to be invoked if the 



lifetime of mass-loss of £ Pup has removed enough of the 
star's envelope to expose nuclear processed material. 

There is some uncertainty r egarding the distanc e to 
£ Pup. The spectroscopi c parallax ( Markova et al ]|2004T ) and 
trigonometric parallax l|Perrvman et al.lll997l ) are in good 
agreement (~ 460 pc and 429 j^ pc, respectively). But it 
has also been suggested that £ Pup lies farther away, at 
d ~ 730 pc, where its space motion and age are consis- 
tent with an origin in the Vela Molecular Ridge (|S ahull 19*921 ; 
ISahu fc Blaauwlll993l ). This larger distance implies a larger 
radius and an Ha mass-loss rate that is larger by a factor 
of two. On the other hand, the Hipparcos data has recently 
been r eanalyzed and a sm aller distance - 332 pc - has been 
found (Ivan Leeuwedl2007T ). 

We stress that the adopted distance does not affect the 
X-ray line profile fitting results directly. But it does affect 
the mass-loss rate we derive from our fits via the depen- 
dence of M on R t , and it affects the fiducial mass-loss rate 
to which we compare the value we derive from the X-ray pro- 
files in this paper. The Ha and radio mass-loss rates scale as 
M oc d 1 ' 8 and the mass-loss rate we derive from the profile 
fitting results scales as M oc d, so the ratio scales only as the 
square root of the distance. Thus, any change in the distance 
will not strongly affect the discrepancy we find between the 
fiducial mass-loss rate and the one we derive from the X- 
ray line profiles. The radius we use for our mass-loss rate 
calculation in this paper assumes the spectroscopic parallax 
distance of 460 pc, which is also assumed for the fiducial Ha 
mass-loss rate determination. 

Detailed spectral synthesis has been carried out from 
the UV to the IR to determine the stellar and wind prop- 
erties of £ Pup, which we list in Table [1] Most of these pa- 
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Table 1. Stell ar and wind parameters adopted from 
IPuls et al.l l|2006h 



parameter 



value 



Distance 


460 pc 


Mass^ 


53.9 Mq 


T e g 


39000 K 


R* 


18.6 R 


frotsini^ 


230 km s" 1 


Voo 


2250 km s" 1 





0.9 




8.3 X 10" 6 M yr" 1 




4.2 X 10" 6 M yr" 1 



a From lRepolust et al.l j2004T) . 

b From lGlebocki et al.l feOOOT). 

c Unclumped value from lPuls et al.l l|2006h . 
^ Also from IPuls et"aT 

I 1120061) . but the minimum clumping 
model, in which the far wind, where the radio emission arises, 
is unclumped, but the inner wind, wher e the Hq is prod uced is 
clumped. Note that the methodology of 1 Puis et al .1 l i2006t) only 
enables a determination to be made of the relative clumping 
in different regions of the wind. This value of the mass-loss 
rate, therefore, represents an upper limit. 



rameters are taken from IPuls et all l|200d ). There is a ranee 
of wind property determinations in the extensive literature 
on £ Pup. The t erminal velocity of the win d may be as low 
as 2200 km s " 1 dLamers fe Leithererlll993h . and as high as 
2485 km s _1 (|Prinia et alJll990h . though we adop t the de- 
termination by the Munich group l|Puls et alJ^OOrj ). of 2250 
km s~\ as our standard. 

Mass-loss rate determinations vary as well. This is 
partly because of the uncertainty in the distance to £ Pup. 
But, it is also the case that each mass-loss rate diagnostic is 
subject to uncertainty: density-squared diagnostics like Hq 
and free-free emission are affected by clumping, no matter 
the size scale or optical depth of the clumps. Mass-loss rates 
from UV absorption lines are subject to uncertain ionization 
corrections. In the last few years there have been attempts 
to account for clumping when deriving mass-loss rates from 
both density-squared diagnostics and UV absorption diag- 
nostics. We list two recent Hq mass- loss rate determina- 
tions in the table, one that assumes a smooth wind and one 
that parameterizes small-scale clumping using a filling fac- 
tor approach. The X-ray line profile diagnostics of mass-loss 
rate that we employ in this paper are not directly affected 
by clumping; although very large scale porosity (associated 
with optically thick clumps) can affect the profiles, as we 
have already discussed. 

The star shows period ic variability in various UV wind 
lines |Howarth et al.lTl995l ) as well as Ha |Berghoefer et ail 
1996). Its broad-band X-ray properties are nor mal for an O 
star, with L x « 10 _7 I/Boi and a soft spectrum (|Hillier et al.l 
1 19931 ). dominated by optically thin thermal line and free- free 
emission from plasma with a temperature of a few million 
degrees. The emission measure filling factor of the wind is 
small, roughly one part in 10 3 . Weak soft X-ray variability, 
with an amplitude of 6 p ercent, and a period of 16.7 hr, was 
detected with ROSAT (|Berghoefer et alJll996h . This low- 
level variability appears not to affect the Chandra data. 



4 EMISSION LINE PROFILE MODEL 
FITTING 

4.1 The Model 

The X-ray emission lin e profile model w e fit to each line 
was first de scribed by [ Owocki fc Cohenl d200ll). b uildin g 
on work by iMacFarlane et al.l (|l99ll ) and llgnacd (|200 if ). 
It is a simple, spherically symmetric model that assumes 
the local emission scales as the ambient density squared 
and that the many sites of hot, X-ray emitting plasma are 
smoothly distributed throughout the wind above some on- 
set radius, R a , which is expected to be several tenths of 
a stellar radius a bove the photosphere in the line-driven in- 
stability scenari o (Owock i et al. 1988; lFeldmeier et al. I ll997l ; 
iRunacres fc Owockil 120021 ). Attenuation of the emitted X- 
rays occurs in the bulk, cool (T « T e fj) wind component 
via photoelectric absorption, mainly out of the inner shell 
of elements N through Si and also out of the L-shell (n = 2) 
of Fe. Singly ionized helium can also make a contribution 
at long wavelengths. We assume that the atomic opacity of 
the cool wind, while a function of wavelength, does not vary 
significantly with radius. This is confirmed by our non-LTE 
wind ionization modelling, discussed in §5.1. We further as- 
sume a beta- velocity law, v — Vao(l — R*/r)^, for both wind 
component s, with Vqq = 22 50 km s _1 as given by UV ob- 
servations (|Puls et al.ll2006l ). The local velocity controls the 
wavelength dependence of the emissivity, the local optical 
depth governs the wavelength-dependent attenuation, and 
the density affects the overall level of emission. The first 
two of these effects can be visualized in Fig. [2] 

We cast the expression for the line profile first in spher- 
ical coordinates, but evaluate some of the quantities explic- 
itly in terms of ray coordinates, with the origin at the center 
of the star and the observer at z = oo. We integrate the spe- 
cific intensity along rays of given impact parameter, p, and 
then integrate over rays. Integrating over the volume of the 
wind, we have: 



, 2 —T<u,r) 

r)r e 



dr. 



(1) 



where L\ is the luminosity per unit wavelength - it is the 
X-ray line profile. The angular coordinate fi = cos 6, and 
r) is the wavelength-dependent emissivity that accounts for 
the Doppler shift of the emitting parcel of wind material 
(which is completely determined, under the assumptions of 
spherical symmetry and the velocity law, according to its 
location, (fi, r) or (p, 2)). The emissivity has an additional 
radial dependence due to the fact that it is proportional 
to the square of the ambient plasma density. The optical 
depth, r, is computed along a ray, z = jj,r, for each value of 
the impact parameter, p = y/l — fi 2 r, as 



r(fi,r) = t(p,z) = 



K,p(r')dz', 



(2) 



w here the dummy radial coordinate is given by r' = 
\J z' 2 + p' 2 . The opacity, k, does not vary significantly across 
a line (recall it is due to continuum processes - the strong 
wavelength dependence across a line profile arises purely 
from the geometry indicated in Fig. [2} . Using the continuity 
equation and the beta-velocity law of the wind, we have: 
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Figure 2. A visualization of the wind Doppler shift and op- 
tical depth - two effects that govern the observed, broadened 
and asymmetric line shapes. The observer is on the left, and the 
light contours represent the line-of-sight velocity in increments of 
0.2i;oo , with the blue shifts arising in the left hemisphere and the 
red shifts in the right. The star is the gray circle at the center, 
and the inner radius of the wind X-ray emission, R a , is indicated 
at 1.5 R* by the solid black circle. The solid heavy contour rep- 
resents the locus of points with optical depth r = 0.33, and the 
dashed and dotted contours represent r = 1 and 3, respectively. 
The model parameters visualized here are nearly identical to those 
of the best-fitting model for the Ne X Lyct line shown in Fig. [8]- 
Ro = 1.5; T* = 2. 



t(p,z) 



R,dz' 



r' 2 (l - R„/r'y 



(3) 



We account for occupation of the back of the wind by 
the star by setting this optical depth integral to oo when 
p < R* and z < yi?, 2 — p 2 . The constant at the front of 
eq. r* = kM/47tR»v oc , is the fiducial optical depth and 
is equivalent to the optical depth value along the central 
ray, integrated down to the stellar surface, in the case where 
v = ^oo- This quantity, t„, is the key parameter that de- 
scribes the X-ray attenuation and governs the shifted and 
asymmetric form of the line profiles. 

We note that the optical depth integral, while generally 
requiring numerical integration, can be done analytically for 
integer values of /3. We use ft = 1 throughout this paper 
(though we report on tests we did for non-integer f3 values in 
i]4.3|l , and for that value of the parameter, the optical depth 
integral along a ray with impact parameter, p, is given by: 

, _ , 7?*r» / R, 7r 
t(p > R„, z) — ( arctan h — 



z t V 



t(p < R*,z) 



R,T, 

2z« 



log 



Zt 

R,u z 
— arctan arctan — 

zt z t 

R, — z, _R*/i + z, z + z, 
R* + 2» R*fJ, — z* z — z* 



(4) 



(5) 



where z t = \/ p 2 — R* 2 and z t = y R t 2 — p 2 , and the inte- 
gral has been evaluated at z and oo. 

The intrinsic line profile function we assume for the 
emissivity at each location is a delta function that picks 
out the Doppler shift line resonance, 



X] oc 5{X - A [l - fiv(r)/c]}. 



(6) 



This assumption is justified because the intrinsic line width 
is dominated by thermal broadening, which is very small 
compared to the Doppler shift caused by the highly super- 
sonic wind flow. 

Calculating a line profile model, then, amounts to solv- 
ing equations [T] and [3] for a given set of parameters: R a , 
r», the normalization (which determines the overall level of 
rf), and an assumed wind velocity law, described by j3 and 
v oo. This last parameter, Voo, influences the emissivity term 
through its effect on the Doppler shift as a function of radius 
and spherical polar angle. And for our choice of j3 — 1, eqs. 
[4] and [5] replace eq. [3] 

The model produces broad emission lines where the 
overall width (in the sense of the second moment of the pro- 
file) , for an assumed wind velocity law, is governed primarily 
by the parameter R . The closer to the star's surface R is, 
the more emission there is from low-velocity wind material, 
which contributes to the line profile only near line center. 
The value of r* affects the line's blue shift and asymmetry. 
The higher its value, the more blue shifted and asymmetric 
the profile. Large values of r* also reduce the profile width 
by dramatically attenuating the red-shifted emission com- 
ponent of the line. The interplay of the two parameters can 
be seen in fig. 2 of lOwocki fc Cohenl (|200ll ). 

4.2 Fitting the data 

4-2.1 statistical fitting of individual lines 

For each line in the spectrum, our goal is to extract values 
for the two parameters of interest - r» and R - and to place 
formal confidence limits on these values. We begin the anal- 
ysis procedure for each line by fitting the weak continuum 
simultaneously in two regions, one to the blue side of the 
line and one on the red side (but excluding the wavelength 
range of the line itself). We assume the continuum is flat 
over this restricted wavelength region. We then fit the emis- 
sion line over a wavelength range that is no broader than the 
line itself (and sometimes even narrower, due to blends with 
nearby lines, which can induce us to exclude contaminated 
portions of the line in question). The model we fit to each 
line is the sum of the empirical line profile model - described 
by equations [T] [4] and [5] - and the continuum model deter- 
mined from the fit to the two spectral regions near the line. 
Note that the inclusion of the continuum does not introduce 
any new free parameters. The overall model thus has only 
three free parameters: the fiducial optical depth, r* , the min- 
imum radius of X-ray emission, R , and the normalization 
of the line, fn some cases, where lines are blended, we fit 
more than one profile model simultaneously, as we describe 
below, but we generally keep the two main parameters of 
each profile model tied together, and so the only new free 
parameter introduced is an additional line normalization. 

We fit the wind profile plus continuum model to both 
the MEG and HEG data (positive and negative first or- 
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ders) simultaneously, if the HEG data are of good enough 
quality to warrant their inclusion, and to the MEG d ata 
only if they are not. We use the C statistic (|Cas hl ll979h a.s 
the goodness-of-fit statistic. This is the maximum likelihood 
statistic for data with Poisson distributed errors, which these 
photon-counting X-ray spectra are. Note that the maximum 
likelihood statistic for Gaussian distributed data is the well- 
known x 2 statistic, but it is not valid for these data, which 
have many bins with only a few counts, especially in the 
diagnostically powerful wings of the profiles. 

We determine the best-fitting model by minimization of 
the C statistic using the fit task in xspec. Once the best- 
fitting model is found, the uncertainties on each model pa- 
rameter are assessed using the A\ 2 formalism^ outlined in 
chapter 15 of IPress et al.l (|2007T ). which is also valid for AC. 
We test each parameter one at a time, stepping through a 
grid of values and, at each step, refit the data while letting 
the other model parameters be free to vary. The 68 per- 
cent confidence limits determined in this manner are what 
we report as the formal uncertainties in the table of fitting 
results, below. We also examine the confidence regions in 
two-dimensional sub-spaces of the whole parameter space 
in order to look for correlations among the interesting pa- 
rameters. Note that we include an extensive discussion of 
modelling uncertainties in £14.31 

We use the relatively strong and unblended Fe xvn line 
at 15.014 A to demonstrate this fitting process. We show the 
MEG and HEG data for this line, along with the best-fitting 
model (the set of model parameters, r«, R a , and normaliza- 
tion that minimizes the C statistic) in Fig. [3] The best- 
fitting model parameters are: r* = 1.97, R Q = 1.53 R*, and 
a normalization of 5.24 x 10 -4 photons s _1 cm~ 2 . Using the 
AC criterion and testing each of these parameters one at a 
time (while allowing each of the other parameters to vary), 
we find that the 68 percent confidence limits on the fit pa- 
rameters are 1.63 < r* < 2.35, 1.38 < Ro/R* < 1.65, and 
5.04 x 10~ 4 < norm < 5.51 x 10~ 4 . The confidence limits 
should be thought of as probabilistic statements about the 
chance that the true parameter values lies withing the given 
range, given the physical assumptions of the model. 

In Fig. |4] we show 68, 90, and 95 percent confidence lim- 
its in two-dimensional r» , R parameter space. We calculate 
a grid of models (typically 36 by 36), optimizing the other 
free parameters (just the normalization, in this case) at each 
point in the grid, and use values of AC = 2.30,4.61,6.17 
|Press et al.l 120071 ') to define the extent of the confidence 
limits. Plots such as this one are a good means of exam- 
ining correlations between model parameters, in terms of 
their abilities to produce similar features in the line pro- 
files. We can see what the trade offs are between parameters 
in a quantitative way. For example, there is a modest anti- 
correlation between R and r* evident in the figure. Low 
values of R a (shock onset close to the photosphere) reduce 
emission on the line wing relative to the core (because there 
is more emitting material at low velocity). So although low 
values of R (hot plasma as close as 1.15 R*) are allowed 
at the 95 percent confidence limit, they require a large wind 



This criterion is a specific numerical value of AC = Ci — C m i„ 
for model realization i, where C m i n is the C statistic value for the 
best-fitting model. 
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Figure 3. The Fe XVII line at 15.014 A in the MEG (top) and 
HEG (bottom), with the best-fitting model superimposed. We 
have not done any rebinning of the data. The error bars repre- 
sent Poisson, root-N, statistics. The dashed vertical lines indicate 
the laboratory rest wavelength of the emission line, and the two 
dotted vertical lines in each panel indicate the wavelengths as- 
sociated with the Doppler shift due to the stellar wind terminal 
velocity of 2250 km s ■ The model is shown as the thick (red) 
histogram, while the data are shown as (black) solid squares with 
error bars. The fit residuals are shown in the horizontal windows 
below the data, with the same one sigma error bars that are shown 
with the data. 



optical depth, r* «s 3, to compensate. High t» values make 
lines narrower, as small values of 7i do, but they also cause 
lines to be more blue-shifted and asymmetric. So, there is 
some degeneracy between these two parameters, but it can 
be broken for good quality data. We note that the confidence 
limits listed in the table of model fitting results, which are 
for individual parameters considered one at a time, will tend 
to differ somewhat from those inferred from these plots of 
joint confidence limits. 

The value of r* at A = 15 A expecte d from the smooth- 
wind Ha mass-loss rate (|Puls et al.ll2006h is r* = 5.30, using 
the opacity model described in §5.1 (which gives a value of 
k = 37 cm 2 g _1 at 15 A). The best-fitting model with fixed 
t» = 5.30 is shown in Fig. [5] This model does not provide 
a good fit, having AC = 64, implying rejection probabili- 
ties well above 99.99 percent. This is the quantitative basis 
for claims that the X-ray emission lines of O stars in gen- 
eral, and £ Pup in particular, are too symmetric and un- 
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Figure 4. Confidence contours (68, 90, and 95 percent) for the 
model fitting of the the Fe XVII line at 15.014 A. The best fit, 
shown in Fig. Q5] is represented by the filled circle. 
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Figure 5. The Fe XVII line at 15.014 A in the MEG (top) and 
HEG (bottom), with the best-fitting model having r* = 5.30 
superimposed. This is the value implied by the smooth- wind Ha 
mass-loss rate and our wind opacity model. The normalization 
and R were the adjustable parameters of this fit. Even this best- 
fitting model is statistically unacceptable. 




6.60 



6.65 6.70 

Wavelength (A) 



6.75 



6.S0 




0.04 
0.00 
-0.04 



6.60 



6.65 6.70 6.75 

Wavelength (A) 



6.80 



2.0 



1.8 



1.6 



1.4 



1.2 




Figure 6. The MEG (top) and HEG (middle) measurements of 
the Si XIII helium-like complex near 6.7 A, along with the best- 
fitting model. This line complex shows a relatively small degree 
of blue shift and asymmetry, indicative of a low r» value, as is 
expected at short wavelengths, where the wind opacity is smaller. 
Note that there is a separate set of vertical lines — denoting the 
rest wavelength and the Dopplcr shifts associated with the wind 
terminal velocity — for each of the three components of the line 
complex (resonance, intercombination, and forbidden lines, from 
short to long wavelength). In this and the following three figures, 
we also show the 68, 90, and 95 percent confidence limits in r* , 
Ro parameter space (bottom). 
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Figure 7. The derived value of t» is also low for the Mg XII Lya 
line at 8.421 A shown here, but it is modestly higher than the 
shorter wavelength Si XIII complex shown in the previous figure. 
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Figure 8. The Ne X Lya line at 12.134 A shows an intermediate 
degree of blue shift and asymmetry, indicative of an intermediate 
t* value, as is expected at its wavelength, where the wind opacity 
is larger than at the wavelength of the Mg XII Lya line, but not 
as large as at longer wavelengths. Part of the red wing of this line 
has been excluded from the fitting because of a possible blend 
with an iron line. 



shifted to be explained by the standard wind-shock scenario 
dKahn et al.ll200ll;ICassinelli et alj|200ll ; iKramer et aT1l2003l ; 
Oski nova et al.ll2006t l. However, the primary goal of this pa- 
per is to quantify the mass-loss rate by modelling the wind 
opacity and the effects of wind attenuation on all the line 
profiles together. To enable us to do this, we repeat the fit- 
ting procedure described above for all 21 of the lines and line 
complexes in the spectrum that have more than 50 counts. 



4-2.2 fitting helium-like line complexes 

For the helium-like complexes - Ovu, Ne IX, Mg XI, Si 
XIII, and S xv - we fit a modified version of the wind profile 
model in xspec that simultaneously fits three separate pro- 
files with the basic parameters (r* and R a ) tied together. It 
accounts for the altered forbidden-to-intercombination line 
strength ratios due to the effects of photoexcitation out of 
the 2 3 S state, which is the upper level of the forbidden line. 
This model, which was first described in iLeutenegger et al.l 
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Figure 9. The O VIII Lya line at 18.969 A shows a relatively 
large degree of blue shift and asymmetry, indicative of a higher 
T * value, as is expected at longer wavelengths, where the wind 
opacity is larger. We did not include the very weak HEG data in 
the analysis of this line. 

(2006), assumes a spatial distribution of X-ray emitting 
plasma, just as the basic wind profile model does, but al- 
ters the radius-dependent line ratio according to the ultra- 
violet mean intensity computed from an assumed model at- 
mospher^fl This model thus self-consistently accounts for 
the effects of the radial dependence of the individual line 
emissivities on both the line ratio and the profile shapes. 
Although the components of these complexes are blended, 
we can extract useful model parameters and confidence lim- 
its on those parameters by fitting each complex as a single 
entity. 

4-2.3 line blends 

We handle other line blends in a manner similar to the 
helium-like complexes, simultaneously fitting profile mod- 
els with parameters tied together. However some blends - 
composed of lines from different ionization states or different 
elements - are more problematic, as their relative strengths 
are generally more uncertain. In some cases, the blending 
is mild - through a combination of the second line being 

2 TLUSTY O star model jLanz fe Hubenvl 120031 ) with T eff = 
40000 K and logg interpolated between 3.50 and 3.75. 



weak and the overlap region being small - and we can fit 
the stronger of the components reliably by simply excluding 
some of the data. This was the case for the Ne x Lya line 
at 12.134 A, where the extreme red wing is mildly blended 
with a weak iron line. In other cases, like the Fe xvn lines 
at 17.051 and 17.096 A, where the relative intensities of the 
components are well constrained by atomic physics, we ob- 
tain reliable results. But this situation is generally only true 
when the overlap between lines is modest and, especially, 
when both lines arise from the same ionization state of the 
same element, as is the case for these Fe xvn lines. 

There are, however, several blends where the modelling 
of the relative line strengths is simply too uncertain to draw 
any reliable conclusions. A good example of this is the N 
vn Lya line at 24.781 A, which is blended with the N VI 
He/3 line at 24.890 A. For this line blend we fit a series of 
models with two components - one for the Lya line and one 
for the He/3 line - trying different values of their relative 
normalizations, all within a plausible range (of 0.1 to 0.4) 
as impli ed by the Astroph ysical Plasma Emission Database 
(aped) l|Smith et alJl200lf ). We found values for the fiducial 
optical depth, r», ranging from less than 1 to more than 
4. These constraints are nearly meaningless, and thus we 
exclude the results for this blended line complex from the 
overall analysis described in §5. 

The five line blends that could not be reliably fit are 
indicated in Fig. [1] by the dashed vertical lines between the 
panels. We stress that we fit all five of these complexes with 
multiple profile models, and in each case found that it was 
impossible to put reliable constraints on r* and R , given 
the wide range of possible relative line normalizations. In 
addition to the N vn and N VI blend near 24.8 A, the com- 
plexes that we had to reject include the helium-like neon 
complex near 13.5 A, which is blended with at least seven 
iron lines that have relative strengths which are temperature 
dependent, and the O Ly/3 line at 16.006 A, which is blended 
with four Fe XVIII lines, and two complexes near 11.0 A and 
15.26 A that contain numerous weak lines arising from Fe 
xvi, Fe xvn, and several higher ionization states of iron. 

4-2.4 results 

After eliminating the five line complexes too blended to give 
meaningful results, we are left with sixteen lines and line 
complexes that could be fit with the wind profile model as 
described in the previous subsection and as demonstrated 
on the Fe xvn line at 15.014 A. The results of these fits are 
summarized in Table [2] And we show four more represen- 
tative line fits - spanning a wide range of wavelengths and 
derived values of t* - in Figs. [6] [7] [HI and[9] Note the progres- 
sion in these profiles from fiducial optical depths, r» , close 
to zero at the shortest wavelengths to significantly larger 
values (up to t, = 3) at the longest wavelengths. We sum- 
marize the sixteen derived r* and R values, along with their 
confidence limits, in Fig. 1101 

4.3 Sensitivity of fitting results to modelling 
assumptions 

We have made various assumptions and choices in carrying 
out the line-profile modelling described in the previous sub- 
section. And we therefore have investigated many of these, 
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Table 2. Wind profile model fit results 



ion 


wavelength 
(A) 


T* 




Ro 
(R*) 


h 

normalization" 
(10~ 5 ph cm -2 s _1 ) 


S xv 


5.0387, 5.0648, 5.1015 


0.01+ ■** 

— .01 


1. 




2.56+-^ 

— .36 


Si xiv 


6.1822 


0.49+-" 

— .35 


1 




0.77+" 

— . 14 


Si xiii 


6.6479, 6.6866, 6.7403 


0.42+-" 

— .Id 


1 


50 +: S 

— .04 


11.2+1 


Mg xi 


7.8503 


0.65±|| 


1 


.33+!? 

— .13 


1.33+'" 

— .13 


Mg xii 


8.4210 


1.22+-g 


1 


.34+-J? 


2.951- 


Mg xi 


9.1687, 9.2297, 9.3143 


0.92+ : ^ 


1 


cz+'-Oe 
oo -.06 


17.8+/® 


Ne x 


9.7082 


0.62l ± 5 ° 5 


1 


•™-.19 


0.95±;i» 


Ne x 


10.2388 


1.95f;g 


1 




2.99+2g 


Ne ix 


11.5440 


0.83±-|| 


2 




5.orf 


Ne x 


12.1339 


2.03±i« 


1 


17+' 11 
•*'-.io 


26.91 1 ,. 1 


Fe xvn 


15.014 


IM±™ 


1 


55+ 13 
°°-.12 


52.412-5 


Fe xvn 


16.780 


2.86+-^ 


1 


01 + .61 
ui -.00 


23.11}- 


Fe xvii^ 


17.051, 17.096 


2 .52t;™ 


1 


'^'-.46 


32.7±;-; 


O viii 


18.969 


3 - 02 -'.57 


1 


.I8t£ 


37.011* 


N vii 


20.9099 


4 - 26 -l'.71 


1 


00 + -87 
.»«_ 87 


14.8l?;| 


O vii 


21.602, 21.804 


1.62l 1 7 ^ 3 


2 


CO + .85 
• OJ -.50 


59.9+« 



a Closely spaced doublets in the Lyman series lines and He-like in- 
tercombination lines are fit with a single profile model centered at the 
cmissivity- weighted wavelength of the two components. 
" For the blended lines fit simultaneously, including the He-like com- 
plexes, the total normalization of all the lines in the complex is indi- 
cated. 

c We fit these two blended lines simultaneously, with a fixed normal- 
ization ratio of 0.9. Both line profile components were forced to have 
the same t* and R values. Allowing the intensity ratio to vary be- 
tween 0.8 and 1.0 hardly changed the parameter confidence limits at 
all. 



again using the Fe xvn line at 15.014 A as a test case. In 
this subsection, we report on the sensitivity of our results 
to the following assumptions and choices: background sub- 
traction; determination of the continuum level; exclusion of 
portions of the line due to possible blending; inclusion of 
the weak HEG data; the adopted values of /3 and Voo for the 
wind; and whether to allow the X-ray volume filling factor 
to vary with radius (as parameterized by q in fx oc r~ q , 
wher e the filling factor , fx, c ontributes to the emissivity, r\ 
- see lOwocki fe Cohenl (|200ll 'l1. We will very briefly describe 
those factors that we found to be unimportant, and discuss 
in more detail those that did make a difference. The baseline 
model fitting we describe here is the modelling described in 
the previous subsection for the 15.014 A line, except that we 
fit only the MEG data (so that we may evaluate the effect 
of including the HEG data). 

We examined the default background spectra, which 
were very weak, and also experimented with fitting the 
15.014 A line with and without the background spectrum 
subtracted and found almost no difference in the fit qual- 
ity or fit parameters. We therefore opt to neglect the back- 
ground when fitting each of the lines in the spectrum. The 
sensitivity to the continuum fit is a little greater, but still 
nearly negligible. When we changed the continuum level by a 
factor of two - which is larger than the formal uncertainty on 
the continuum level - none of the parameter values changed 
by more than ten percent. Some lines in the spectrum are 



blended with weaker lines. The cleanest way to handle this 
situation is to exclude the contaminated bins from the mod- 
elling. To test the effects of this, we eliminated 0.03 A from 
the red wing of the 15.014 A line and refit the data. We then 
repeated this experiment eliminating 0.07 A - leaving only 
about two-thirds of the data. Even in this second, extreme 
case, the fit parameters varied by less than ten percent and 
the confidence regions only expanded slightly. 

For most lines, the HEG data is significantly weaker 
than the MEG data. We find for the 15.014 A line that 
including the HEG data changes the best-fitting model pa- 
rameters by, at most, a few percent, but it does tighten the 
confidence limits somewhat. The effect of including the HEG 
data is more significant for the shorter wavelength lines, 
where the effective area of the HEG is larger relative to the 
MEG. There is very little penalty for including the HEG 
data, so we do so for all lines shortward of 16 A. We also 
fit the MEG and HEG data separately for the 15.014 line 
to verify that there are not systematic effects between these 
two spectra; and there are not. The separate fits give re- 
sults that are very similar to each other, with significantly 
overlapping 68 percent confidence limits for all parameters. 

The original lOwocki fe Cohenl l|200ll ) line profile model 
allows for a radially varying filling factor of X-ray emitting 
plasma, parameterized as a power law function of radius. 
Values of the power-law index, q, that differ significantly 
from zero (no radial variation) can cause changes in the line 
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Figure 10. Values of r* (top) and i? (bottom) derived from the 
model fits, shown with their 68 percent confidence limits. Line 
complexes and blends that were fit with multiple model compo- 
nents are represented by only one point. 



profiles that are not insignificant, effectively weighting the 
emission from parts of the wind according to their velocity 
(via the beta-velocity law relationship between velocity and 
radius). However, we find that when we allow q to be a free 
parameter the best-fitting value is generally very close to 
zero. For the representative 15.014 A line, it is q = —0.09, 
and q — is included in the 68 percent confidence range. 
The general result is consistent with that found for this and 
other stars (jKramer et alj|2003l ; ICohen et al1l2006t ). Thus, 
to keep the number of free parameters manageable, we fix 
9 = 0. 

The factors discussed above have a very minor influence 
on the results of the line fitting. However, the remaining 
factors can have a significant effect. 

The velocity-law exponent, /3, affects line profiles for 
two reasons: (1) the velocity law affects the mapping be- 
tween radius and Doppler shifted wavelength, and so affects 
the emission profile; and (2) via the continuity equation, 
it affects the density and so affects the level of both the 
emission and the absorption. Indeed, for our representative 
emission line, when we change the value of ft from 1 to 0.8, 
both T t and R change by 10 to 20 percent. The determi- 
nations of P for £ Pup vary from at least 0.9 to 1.15, and 
so using a value of /? = 1 seems reasonable, especially as it 



speeds the calculation of the line profile model by allowing 
the optical depth integral to be done analytically, so we use 
that value for all the model fitting results reported here. If, 
in the future, a new and more accurate determination of ft 
is made, and it differs significantly from /3 = 1, then the re- 
sults reported in this paper can be scaled according!^. We 
also note that the X-ray emitting plasma and the bulk wind 
that attenuates the X-rays may not necessarily be described 
by the same beta velocity law. However, there is no indepen- 
dent evidence for this, and with the short post-shock cooling 
lengths expected in the relatively dense wind of £ Pup, the 
X-ray emitting plasma in the wind is more likely to have a 
velocity close to the ambient wind velocitjQ. And further- 
more, the observed X-ray emission line widths in £ Pup and 
other early O supergiants are completely consistent with the 
P and Woo values inferred from UV and optical spectroscopy 
of these stars. 

The terminal velocity of C, Pup is relatively well es- 
tablished, with reasonable estimates from several different 
groups that vary by about ±10 percent about our adopted 
value of 2250 km s _1 . However, when we explored the effect 
of varying the terminal velocity in our fitting of wind pro- 
file models to the 15.014 A line, we found that the value of 
t* was quite sensitive to the assumed wind terminal veloc- 
ity, even within this relatively narrow range. This is because 
the blue shift of the line centroid in the dimensionless, scaled 
wavelength parameter, x = (X/X — l)c/voo, depends directly 
on the degree of wind absorption. The same observed pro- 
file appears more blue shifted in scaled wavelength units if 
the terminal velocity is (assumed to be) smaller. Our tests 
with the 15.014 A line show that the best-fitting value for r» 
ranges from 2.16 to 1.35 when we use terminal velocities be- 
tween 2200 km s _1 and 2485 km s _1 . This variation is larger 
than that caused by every other parameter uncertainty and 
assumption we have explored. Thus, while we consider the 
value of Woo = 2250 km s _1 to be quite reliable, future re- 
assessments of this parameter will necessitate a rescaling of 
the optical depth - and mass-loss rate - results we report in 
this paper. 

As a final test, we can treat the terminal velocity as 
a free parameter of the model. This enables us to see what 
value of the terminal velocity is preferred by the X-ray spec- 
tral data themselves. In general, the constraints on Woo , while 
letting the other model parameters be free to vary, were not 
strong. But for the highest signal-to-noise lines in the spec- 
trum, relatively tight constraints could be derived. We show 
the results for fitting the five most useful lines in Fig. 1111 
As the figure shows, these lines are all consistent with our 
adopted value of Voo = 2250 km s~ . This, of course, gives 
us added confidence that the value we use for the model fit- 
ting is reasonable. And, in fact, the small error bars on most 



3 Lowering f} from 1 to 0.8 causes the best-fitting optical depth 
of the Fe XVII line at 15.014 A to go from t* = 1.98 to t* = 1.66. 
If the value of were to be revised upward by a similar amount, 
the values we derive for r» from the line profile fitting would have 
to be revised upward by about 15 percent. The quality of the fits 
with the different values of /3 do not differ significantly. 

4 X-ray emitting plasma is too highly ionized to be effectively 
driven by the photospheric UV radiation field. However, for small 
enough parcels, the ram pressure of the surrounding wind should 
keep the post-shock, hot plasma moving at the ambient velocity. 
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Figure 11. Values of the terminal velocity derived from fitting 
five strong lines with a wind profile model for which Voo was 
allowed to be a free parameter (the other parameters - r*, Ro, 
and the normalization - were allowed to vary as well). The bulk 
wind terminal velocity adopted from the analysis of UV profiles 
is indicated by the solid horizontal line. The cross-hatched area 
represents the 68 percent confidence region for the value of the 
terminal velocity derived from fitting these five points. 



of these determinations also show that significantly smaller 
and larger values are ruled out. The kinematics of the hot, 
X-ray emitting plasma seem to be the same as that of the 
bulk wind for £ Pup. 



5 DISCUSSION 

The most obvious new and significant result of the profile 
model fitting is the wavelength trend in the derived values 
of the fiducial optical depth, r«, shown in the top panel of 
Fig. [10] The value of this parameter, which is proportional 
to both the mass-loss rate and the opacity of the bulk wind, 
increases with wavelength, which is exactly what is expected 
from the form of the atomic opacity. The null hypothesis of 
a constant value of r* is rejected with greater than 99.9 
percent confidence (xt — 5.4 for 15 degrees of freedom). We 
therefore fit a model of wavelength-dependent r» , in which 
the wavelength dependence derives entirely from the atomic 
opacity, k(X). 

While it may seem obvious that there should be a trend 
in the fiducial optical depth with wavelength, this result is 
quite significant, in that a presumed lack of such a trend is 
the basis for claims that large-scale clumping and the as- 
sociated wind porosity are the cause of t he smaller than 
expe cted profile blue shifts and asymmetry (|Oskinova et al.l 
120061 ). In the following subsections, we show how a realistic 
wind opacity model naturally explains the observed wave- 
length trend, and then how such a model can be used to 
make a quantitative determination of the mass-loss rate of 
CPup. 



5.1 The opacity model and the mass-loss rate 
determination 

The opacity model depends on the abundances and, to 
a lesser extent, the ionization balance of the bulk stellar 



wind (i.e. the cooler, unshocked component). The dominant 
source of opacity is photoelectric absorption from the K- 
shell of abundant elements between N and Si, and also the 
L-shell of Fe. We have computed a wind opacity m odel using 
Cmfgen jrlillier fe Milleijll998l ; [Zsarg6 et al.ll2010l ) to model 
the ionization balance an d using atomic cross sections from 
IVerner fe Yakovlevl j 19951 ). The model is constrained by UV 
and optical spectra, so the abundances are de rived directly 
from observat ions. Details are p rovided in iBouret et al.l 
|2010l ) and in IBouret et all ((2008) where the overall mod- 
eling is briefly described and excellent fits to Ha and P v 
profiles are shown. Specifically, it is found that Yh c = 0.16 
((Z/Z0)hc = 1.88 expressed as a fraction of the solar abun- 
dance), (Z/Z Q )c = 0.08, (Z/Z ) N = 5.0, (Z/Z Q )o = 0.20, 
and (Z/Z0)fc — 1.0, where the reference sol a r abu ndances 
are taken from lAsplund. Grevesse. fe Sauvall (|2005l ). These 
abundances are consistent with those derived from inde- 
pendent analysi s by the Munich group (J. Puis, private 
communication; IPauldracbl (|2003l )). Additionally, the low 
oxygen abundance is consistent with the v alue found from 
model ling the X-ray spectrum (0.30 ± 0.43, Zhe kov fe Pallal 
(2007)). These authors also find a high nitrogen abundance 
of 3.2 ± 0.6, only slightly lower than the value we adopt 
here. Note that we ha ve scaled the abundances reported by 
Zhekov fe Pallal d2007Tl to the reference solar abundances of 



Asplund. Grevesse. fe Sauva 3 |2005l ) 



We show the wind opacity model, using our adopted 
abundances and the ionization balance from the cmfgen 
modelling, at a sing le radius (r = 1.8 i?*J3 in Fig. 1121 
along with a solar-abundance model. The opacity is lower 
at most wavelengths in the cmfgen model primarily be- 
cause the total abundance of metals (and most crucially the 
su m of carbon, nitrogen, and oxygen) is subsolar (0.53 of 
the I Asplund. Grevesse. fe Sauvall |2005) value). We refer to 
this opacity model, based on the cmfgen modelling and the 
abundances derived from the UV and optical spectra, as the 
"sub-solar metallicity" model in the remainder of the paper. 

Using either of these models of the opacity, and values 
for the stellar radius and wind terminal velocity from Table 
[JJ we can construct a wavelength-dependent model of r», 
for which the mass-loss rate is the only free parameter. Fits 
with both the sub-solar metallicity wind opacity model and 
the solar abundance model are good (xt ~ 0.6 for the sub- 
solar metallicity model and xt ~ 0.8 for the solar abundance 
model), although a higher mass- loss rate of M = 3.50 x 10 -6 
Mq yr _1 is found with the sub-solar metallicity model, due 
to its lower overall opacity. The solar abundance opacity 
model, which should provide a lower limiting case, gives M = 
1.90 x 10 -6 Mq yr _1 . The formal uncertainties on these 
derived mass-loss rates, due solely to the finite error bars on 
the individual t* determinations, are about 10 percent. 

The best-fitting r* model, using the sub-solar metallic- 



5 We note that there is very little variation in the opacity with 
radius between 1.1 R* and roughly 4 R* (at least at wavelengths 
where we analyze lines, below the nitrogen K-shell edge near 26 
A). By 5 -R* the overall opacity is about twenty percent higher, 
and by 11 ij* it is about a factor of two higher. The increasing 
opacity with radius is due to the larger fraction of singly ionized 
helium in the outer wind. But the wind density is so low at these 
distances that the outer wind does not contribute significantly to 
the X-ray optical depth. 
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Figure 12. The wavelength dependent opacity of the wind of 
£ Pup computed using CMFGEN modelling of the ionization bal- 
ance and sub-solar abundances derived from UV and optical spec- 
tra (solid), along with a solar-abundance opacity model (dotted). 
Note the prominent K-shell edge of oxygen near 20 A in the solar 
abundance model. In the sub-solar metallicity opacity model, this 
decrement is much more modest, due to the underabundance of O 
and overabundance of N. The overall reduction in the opacity at 
most wavelengths in the CMFGEN model is the result of its overall 
subsolar metallicity and not its altered CNO abundances. 

ity opacities and the best-fitting mass-loss rate, is shown in 
Fig. 1131 along with the r, model computed using the smooth- 
wind Ha mass-loss rate, M = 8.3 x 1CT 6 Mq yr" 1 . The 
best-fitting mass-loss rate is almost a factor of three loweijf]. 
If solar abundances are assumed for the opacities, the fac- 
tor is more than four. The best-fitting versions of these two 
models are compared in Fig. 1141 and have a very similar 
shape, implying that even with better quality Chandra data 
it would be difficult to distinguish them based on the X-ray 
data alone. We stress, though, that the abundances of £ Pup 
are certainly not solar. We present this model only for com- 
parison with the sub-solar metallicity opacity model, and as 
a limiting high opacity case. 

Taking a closer look at the atomic opacity, we can see 
in the preceding three figures that the most leverage regard- 
ing the wavelength dependence of the opacity, and hence of 
t* , comes at the shortest wavelengths, below the Ne K-shell 
edges near 13 A. The Fe and Ne edges and the low O abun- 
dance conspire to make the opacity rather flatter than the 
generally expected k oc A 3 relationship seen from individual 
elements' photoionization cross sections. Most of the strong 
lines in the MEG spectra of O stars are between 12 and 18 A, 
where the opacity is relatively constant. This demonstrates 

6 The mass-loss rate we derive here from the X-ray line profiles 
is nearly identical to the Ha a nd radio mass-loss rate determined 
bv lLamers fe Leithererl jl993T ). al though this is purely coinciden- 
tal. Several systematic errors in lLamers fc Leithererl il993h re- 
duce their smooth-wind mass-loss rate determination, compared 
to more modern estimates such a s the one we em ploy for the 
fiducial, unclumped mass-loss rate llPuls et al.ll200d). Th e factors 
that lead to the low value in lLamers^TTeTth*ere^ ~ Jl993h include 
assuming that helium is fully ionized in the outer wind, ignor- 
ing departures from LTE and assuming a lower temperature than 
more modern analyses use, and assuming that the Ho line is op- 
tically thin. 
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Figure 13. Values of r* derived from the line-profile model fits, 
shown as points with error bars (same as the top panel of Fig. 
I10I I. The value of r* expected if the mass-loss rate is 8.3 X 10~ 6 
Mq yr — 1 is shown as the upper dashed curve. Treating the mass- 
loss rate as a free parameter, the best fit value of 3.50 X 10~ 6 
Mq yr — 1 is shown as the lower, varying solid curve. This model 
provides a formally good fit. And both of these models of the 
wavelength-dependent t* use the sub-solar metallicity (cmfgen) 
opacity model. The horizontal dash-dot line is the best-fitting 
constant r* model, as would be expected for a highly porous 
wind, where the effective opacity is completely determined by the 
macroscopic, physical cross sections of the optically thick clumps. 
It does not provide a good fit to the data. 
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Figure 14. The best-fitting model, with M = 3.50 X 10~ 6 
Mq yr , shown in Fig. 1131 is shown here again, but this time it 
is compared to the best-fitting solar abundance t* model (dotted 
curve). The fits are of similar quality, while the solar abundance 
model has a lower mass- loss rate (M = 1.90 X 10~ 6 Mq yr" 1 ) to 
compensate for its higher overall opacity. 

the need for the use of realistic wind opacity models when 
interpreting trends in grating spectra of O stars and explains 
why the wavelength trend of r* values was not noted in the 
initial studies. 

Furthermore, the paucity of useful emission lines long- 
ward of the O K-shell edge makes it difficult to discriminate 
among various wind opacity models, although in principle, 
lines longward of this edge could enable us to diagnose the al- 
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tered CNO-processed abundances with some certainty. And 
emission lines longward of the N K-shell edge near 26 A 
would be especially useful, but there are none in the Chan- 
dra spectrum. The N VII Ly/3 line at 20.910 A is quite weak 
and does not provide a strong constraint on r„, although 
it does favor the sub-solar metallicity opacity model. The 
longest wavelength line which we are able to reliably fit is 
the helium-like O vn complex near 21.8 A. We fit the res- 
onance and intercombination lines simultaneously (the for- 
bidden line is not present due to 2 3 S — 2 3 P photoexcitation 
by the photospheric UV field), with the profile parameters 
t* and R tied together for the two lines. However, the res- 
onance line in this complex may b e subject to resonance 
scattering (|Leutenegger et all I2007T ) - it may be optically 
thick to its own radiation (as distinct from the effects of 
continuum opacity of the overlying wind that leads to the 
observed skewness and blue shifts in all of the line profiles). 
Resonance scattering tends to make broadened, asymmet- 
ric, and blue shifted lines more symmetric, and thus the 
t* value we derive from fitting this complex may be un- 
derestimated. If this is the case, then this line complex too 
would favor the sub-solar metallicity wind opacity model, 
as shown in Fig. 1141 We also note that the only other line 
of the sixteen we analyze that is likely to be optically thick 
to resonance scattering is the O vn Lya line at 18.969 A, 
so the r* determination for that line may also be somewhat 
underestimated. 

We also can see from a careful inspection of the opac- 
ity model that the mass-loss rate determination from fitting 
a set of t» values is mostly sensitive to the cross section 
contributions from N, O, and Fe. Alterations of O and N 
abundances due to CNO processing will have only a mod- 
est effect on the results, however. The sum of the contri- 
butions of N and O (as well as He and C) is what affects 
the overall opacity level between about 15 A and 20 A, 
with Fe - and to a lesser extent, Ne - making a signifi- 
cant contribution at shorter wavelengths. This demonstrates 
that accurate determinations of abundances for O stars are 
perhaps the biggest factor in enabling the determination of 
dumping-independent mass-loss rates from high-resolution 
X-ray spectra. But when fitting a large ensemble of lines 
that span a relatively wide range of wavelengths, knowing 
the overall metallicity is probably sufficient, although in- 
cluding a realistic mixture of elements (and thus absorption 
edges) is important too. 

5.2 Sources of Uncertainty in the Mass-Loss Rate 
Determination 

The uncertainty in the mass-loss rate determination we have 
found from the fits to the ensemble of t* values, derived from 
fitting the individual line profiles, come from three sources. 
The first is the formal uncertainty on the mass-loss rate 
model that stems from the uncertainties on the individual 
line profile fits (represented by the error bars on the r* points 
in Fig.[l3l for example). For the sub-solar metallicity opacity 
model the 68 percent confidence limit range on the fitted 
mass-loss rate extends from 3.25 to 3.73 x 10 -6 Mq yr -1 , 
representing an uncertainty of a little less than 10 percent 
on the best-fitting value of 3.50 x 10" 6 Mq yr" 1 . 

The second source of uncertainty arises from our imper- 
fect knowledge of the wind terminal velocity (and, most im- 



portantly, the terminal velocity of the X-ray emitting plasma 
itself). However, as we have shown (see Fig. Hip , the data 
themselves indicate that our adopted terminal velocity of 
i>oo = 2250 km s _1 is well supported. Three of the lines we 
show in that figure have best-fitting terminal velocity val- 
ues near 2350 km s _1 , which is also the terminal velocity 
derived from a careful analysis of the UV line profiles by 
lHaserl (|l995T ). When we refit the representative Fe xvn line 
at 15.014 A using this higher terminal velocity, we found a 
reduction in our derived r* value of 15 percent. If this scal- 
ing holds for all lines, then using this slightly higher value of 
the terminal velocity will lead to a downward revision of our 
derived mass-loss rate of about 15 percent. (Note that the 
terminal velocity enters into the denominator of the expres- 
sion for t*, and that will mitigate this adjustment slightly.) 
Similar considerations pertain to our assumption about the 
wind velocity parameter, /3. 

The third, and largest, source of uncertainty is due 
to the abundances. We estimate that the abundances de- 
rived for £ Pup from the analysis of UV and optical data 
have a precision of only about a factor of two (jBouret et al.l 
l20ld ). However, we note that they are in good agree- 
ment with the indepen dent, X-ray-based determination from 
IZhekov fe Pallal (|2007T ). providing additional confidence as 
to their accuracy. We can see from the comparison of the 
sub-solar metallicity model to the solar abundance model 
that the mass-loss rate varies by about a factor of two be- 
tween these two assumed opacity models, although the solar 
abundance model is included in our analysis not so much as 
a realistic alternate model, but simply as a plausible up- 
per bound to the atomic opacity; the sub-solar metallicity 
model is more realistic due to the constraints on it provided 
by observa tions in other wavelength bands (|Bouret et al.l 
I2008l . l2010h , and of course, the evolved nature of £ Pup im- 
plies that we should not expect to find solar abundances in 
its wind. However, the overall CNO metallicity of 0.53 so- 
lar is lower than expected, giv en the uniformly sola r abun- 
dances in nearby massive stars (|Przvbilla et al.|[2008T ). Thus, 
a conservative estimate of the allowed range of the mass-loss 
rate of ( Pup derived from the X-ray line profile fitting is 
roughly 2 to 4 x 10 -6 Mq yr" 1 , with our best estimate be- 
ing 3.50 x 10 -6 Mq yr~ . This mass-loss rate is only a lit- 
tle lower t han the maximu m mass-loss rate of 4.2 x 10~ 6 
Mq yr -1 (|Puls et al.l I2006T ). implying a small amount of 
clumping in the outer wind, and a small adjust ment to the 
clump ing factor in the inner wind determined bv lPuls et al.l 
|2006h . 

Any future modification to the accepted abundances of 
( Pup could lead to a change in the mass-loss rate implied 
by our X-ray line profile analysis. To a good approximation, 
such a change would simply involve scaling the X-ray line 
profile mass-loss rate by the reciprocal of the change in the 
overall metallicity, for reasons discussed at the end of the 
previous subsection. 

5.3 Location of the X-ray Emitting Plasma 

The analysis of the sixteen lines and line complexes in the 
Chandra spectrum of ( Pup also enables us to derive values 
of the onset radius of the wind-shock X-ray emission from 
the profiles. These results are shown in the lower panel of 
Fig. 1101 and are completely consistent with the expectations 
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of the wi nd-shock structure induced by the line-driven in - 
stability l|Feldmeier et al.lll997l ; iRunacres fc Owockill2002l) . 
That is, an onset radius of R ~ 1.5 R* (from a weighted 
fit to the results from the sixteen fitted lines and line com- 
plexes; with an uncertainty of 0.1 i?*). We have searched for 
a trend with wavelength in these values and found non^Z). 
Thus, the simplest interpretation is that there is a universal 
radius of the onset of X-ray emission and it occurs near 1.5 
R* (half a stellar radius ab ove the photosphere) . This result 
had already been noted by IKramer et al l l|2003l ). though we 
show it more robustly here. This sam e result can also b e 
seen in the late O supergiant £ Ori (Co hen et al.l [2006). 
And this result is also consistent with the joint analysis 
of X-ray line profile shapes and helium-like forbidden-to- 
intercombination l i ne rat ios for four O stars as described by 
iLeutenegger et alJ 120061 ). 



5.4 Comparison with Previous Analyses 

Finally, let us consider why we have found a trend in wave- 
length for the fiducial optical de pth values, t*, derive d from 
the same Chandra data that led lKramer et all (|2003) to re- 
port that there was no obvious trend. The two biggest fac- 
tors leading to this new result are our more careful assess- 
ment of line blends, and our inclusio n of many weak, bu t 
important, lines at short wavelength. iKramer et all <l2003h 
included only one line shortward of the Ne x Lya line at 
12.134 A, whereas we report on nine lines or line com- 
plexes in_Uns_j^nge ( inclu ding two helium-like complexes, 
which IKramer et all (|2003r i excluded from their analysis). 
While many of these lines are weak and do not provide very 
strong constraints when considered individually, taken to- 
gether, they are high l y stat istically significant. As for line 
blends. IKramer et al.l (|2003T ) included the N vn Lya line at 
24.781 A and the Fe xvn complex near 15.26 A, both of 
which we have determined are too blended to allow the ex- 
traction of reliable information about their intrinsic profile 
shapes. Furthermore, we properly account for the blended 
Fe xvn lines at 17.051 and 17 . 096 A , fitting them simul- 
taneously, while IKramer et al.l l|2003h fit them as a single 
line. 

Nonetheless, if we exclude the blended 17.05, 17.10 
A lines, our r« va lues for each line also analyzed by 
IKramer et all (|2003h are in agreement to within the error 
bars. Similarly, fo r five unblended l ines i n the analysis of 
the same data by lYamamoto et all (|2007T ) , we find consis- 
tent results. In fact, the wavelength tre nd of r« is fully 
consis tent with the r* values found by lYamamoto et all 
(120071 1 . but there were not enough lines in that study for 
the trend to be unambiguously detected. The seven addi- 
tional lines a nd line complexes th a t are not analyzed i n any 
other study iKramer et all 120031 ; ILeutenegger et all 120061 ; 
lYamamoto et al.l 120071 ') but which we analyze here are cru- 
cial for mapping out the wavelength dependence of r*. 

An additional factor that enabled us to determine 



7 An unweighted fit of an assumed linear trend shows a modest 
increase with wavelength, but that result is significant at only the 
one sigma level, and when we perform a weighted fit — with the 
weights inversely proportional to the uncertainties on the individ- 
ual measurements - the significance is less than one sigma. 



that the wavelength trend in r» is consistent with that 
expected from the form of the atomic opacity of the 
wind is our use of a detailed model of the wind opacity. 
It is relatively flat over much of the wavelength range 
encompassing the strong lines in the Chandra spectrum. 
Specifically, from about 12 A to about 18 A, the presence 
of successive ionization edges makes the overall opacity 
roughly flat. Thus, for a trend to be apparent, short 
wavelength lines have to be included in the analysis. 
Previous studies have computed wind X-ray opacities - 
several based on detailed non-LTE wind mo deling - and 
used t hem for the analys is of X -r ay spectra dHillier et al.l 
I l993l; iMacFarlane et"atl Il994l; IPauldrach et al.l 11994 



Cohen et alJ Il996l; iHillier fc Milled Il998l; IWaldron et all 



19981; IPauldrach et all 1200 ll ; lOskinova et"afl 120061 ; 
Krticka fc KubatJ 12009 ) . But the present study is the 



first to demonstrate the importance of the combined effect 
of multiple edges in flattening the opacity in the middle of 
the Chandra bandpass. And it is the first to explore the 
sensitivity of the fitting of a high-resolution X-ray spectrum 
to the assumed wind opacity. 

Finally, the mass-loss rate reduction derived here is 
only a little less than a factor of three, while earlier analy- 
ses suggested that, without porosity, much larger mass-loss 
rate reductions would be required to e xplain the only mod- 
estly shifted and asy mmetric profiles (|Kramer et al.ll2003l ; 
lOskinova et alJ 2006). Here too, the wind opacity is key. The 
overall opacity of the wind is significantly lower than had 
been previously assumed, implying that the mass-loss rate 
reduction is not as great than had been assumed. Again, this 
is primarily due to the significantly sub-solar abundances 
(especially of oxygen) in £ Pup. 



6 CONCLUSIONS 

By quantitatively analyzing all the X-ray line profiles in the 
Chandra spectrum, we have determined a mass-loss rate of 
3.5 x 10~ 6 Mq yr -1 , with a confidence range of 2 to 4x 10 -6 
Mq yr -1 . Within the context of the simple, spherically sym- 
metric wind emission and absorption model we employ, the 
largest uncertainty arises from the abundances used in the 
atomic opacity model. This method of mass-loss rate deter- 
mination from X-ray profiles is a potentially powerful tool 
for addressing the important issue of the actual mass-loss 
rates of O stars. Care must be taken in the profile analysis, 
however, as well as in the interpretation of the trends found 
in the derived r, values. It is especially important to use a 
realistic model of the wind opacity. And for O stars with 
weaker winds, especially, it will be important to verify that 
the X-ray profiles are consistent with the overall paradigm 
of embedded wind shocks. Here, an independent determina- 
tion of the terminal velocity of the X-ray emitting plasma 
by analyzing the widths and profiles of the observed X-ray 
lines themselves will be crucial. In the case of ( Pup, we have 
shown that the X-ray profiles are in fact consistent with the 
same wind kinematics seen in UV absorption line spectra of 
the bulk wind. And the profile analysis also strongly con- 
strains the onset radius of X-ray production to be about 
r = 1.5 R*. 

An additional conclusion from the profile analysis is 
that there is no need to invoke large scale porosity to ex- 
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plain individual line profiles, as the overall wavelength trend 
is completely consistent (within the measurement errors) 
with the wavelength-dependence of the atomic opacity. The 
lower-than-expected wind optical depths are simply due to a 
reduction in the wind mass-loss rate. This modest reduction 
is consistent with other recent determinations that account 
for the effect of small-scale clumping on density-squar ed di- 
agnostics and ionization corrections l|Puls et al.l l2006) . 
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